Elliptic flow in U+U collisions at y/s^N = 200 GeV and in Pb+Pb collisions 
at a/saT/v = 2.76 TeV: Prediction from a hybrid approach 
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We predict the elliptic flow parameter vi in U+U collisions at ^Jsnn = 200 GeV and in Pb+Pb 
collisions at ^/snn = 2.76 TeV using a hybrid model in which the evolution of the quark gluon 
plasma is described by ideal hydrodynamics with a state-of-the-art lattice QCD equation of state, 
and the subsequent hadronic stage by a hadron cascade model. 

PACS numbers: 25.75.-q, 25.75.Nq, 12.38.Mh, 12.38.Qk 



One of the major discoveries at Relativistic Heavy 
Ion Collider (RHIC) in Brookhavcn National Laboratory 
(BNL) was that for the first time in relativistic heavy 
ion collisions, the elliptic flow appeared to be as large 
as an ideal hydrodynamic prediction Since viscosity 
and any other dissipative effects vanish in ideal hydro- 
dynamics, and tiny viscosity requires a strong coupling 
of constituents (quarks and gluons in our case), this dis- 
covery established the new paradigm of strongly coupled 
quark gluon plasma @, Q . 

In noncentral collisions, rcscatterings of the created 
particles convert the initial spatial anisotropyof the reac- 
tion zone to anisotropic particle distribution [J] . Ideal hy- 
drodynamics predicts that the ratio of these anisotropics 
is V2/e ^0.2 almost independent of centrality at the 
RHIC energies Here V2 is the second Fourier coeffi- 
cient of the azimuthal distribution of final particles, and e 
is the initial eccentricity of the produced matter. On the 
other hand, in the dilute regime kinetic theory predicts 
i>2 to be proportional to the particle multiplicity per unit 
rapidity, dN/dy [f| 01 • Thus, the response of the system, 
vijs, provides information about the transport proper- 
ties of the QGP. Experimentally Vi/e is seen to increase 
with increasing transverse density (l/S)dN/dy @, Q, 
where S is the transverse area of the collision zone, until 
it reaches the so-called hydrodynamic limit, t^/e ~ 0.2, 
in central Au+Au collisions at RHIC. With the agree- 
ment of the hydrodynamical prediction of the particle 
mass dependence of V2(pr) with the data [Tl|, [l2T |. 
this is considered as evidence for the discovery of the 
perfect-fluid nature of the QGP 

After observing the increase of vijs with increasing 
transverse density, it is natural to ask what happens if 
the transverse density increases beyond that achieved at 
RHIC [l3| . Will it saturate to the value observed at 
RHIC, as expected if the system behaves like a perfect 



fluid, or will it keep increasing? One suggested way to 
extend the transverse density is to perform uranium- 
uranium collisions [13| |. Since uranium nuclei are de- 



formed and larger than gold nuclei, one can expect large 
transverse density with finite eccentricity in the body- 
body collisions at vanishing impact parameter^]. Some 
Monte Carlo studies show that, even though one cannot 
control the orientation of colliding nuclei, events with 
high multiplicity though finite eccentricity can be se- 
lected in the usual triggering process (T3, Il5l4l8| . An- 
other way to extend the transverse energy is to increase 
the collision energy to generate more particles in colli- 
sions. This is going to happen very soon in the Large 
Hadron Collider (LHC) heavy ion program. In this Rapid 
Communication we predict elliptic flow parameters both 
in U+U collisions at RHIC and Pb+Pb collisions at LHC 
using a hybrid model based on ideal hydrodynamics and 
hadron cascade. 

We describe space-time evolution of the QGP by ideal 
hydrodynamics [19| with the recent lattice QCD equation 
of state [2(J. After expansion and cooling, the system 
turns into hadronic matter. We switch from hydrody- 
namics to a kinetic approach at a switching temperature 
T sw and employ a hadronic cascade model, JAM [2l[, to 
describe the subsequent space-time evolution of hadronic 
matter. 

Our EoS, s95p-vl.l, is a slightly modified version of the 
s95p-vl EoS presented in [20|. It interpolates between 
hadron resonance gas at low temperatures, and recent 
lattice QCD results by the hotQCD collaboration (22l.[23j 
at high temperatures in the same way as s95p-vl , but the 
hadron resonance gas part contains the same hadrons and 
resonances as the JAM hadron cascade [2l[. The details 
of the interpolating procedure are explained in [2fjj and 
the parametrization and EoS tables are available at [24| . 

For initial conditions, we employ two Monte-Carlo 
approaches to simulate collisions of two energetic 
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FIG. 1. (Color online) (a) Initial state eccentricity e par t and (b) V2 as a function of N pa , Tt in Au+Au and U+U collisions at 
^snn = 200 GeV. Experimental data of V2 in Au+Au collisions are obtained by PHOBOS Collaboration [34j |. 



nuclei: Monte-Carlo Glauber (MC-Glauber) model 
[IU and Monte-Carlo Kharzeev-Levin-Nardi (MC-KLN) 
model [26| . In the MC-Glauber model, one calculates the 
number of participants iV pa rt and the number of binary 
collisions N co u for a given nuclear density distribution. 
We model the initial entropy distribution in hydrody- 
namic simulations as a linear combination of the number 
distribution of participants p part = -gaff an d that of 
binary collisions p co \\ = d J^ coH in the transverse plane: 
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We generate the number distributions in an event-by- 
event basis, align them to match the main and sub axes 
of the ellipsoids, and average over many events for a given 
ccntrality bin to obtain a smooth distribution 27 1. The 



eccentricity of the initial profile is then evaluated with 
respect to participant plane, £ part [Hj]. We do the cen- 
trality cuts according to the Af part distribution from the 
MC-Glauber model instead of using the optical Glauber 
limit as was done in Ref. [27| . The free parameters of 
the model, the mixing parameter a = 0.18 and the pro- 
portionality constant in Eq. ((T|) , arc chosen to reproduce 
transverse momentum spectra for pions, kaons, and pro- 
tons from central (0-5%) to peripheral (70-80%) events in 
Au+Au collisions at v / s]v r /v = 200 GeV obtained by the 
PHENIX Collaboration (29[. We also choose the switch- 
ing temperature as T sw = 155 MeV to describe the rela- 
tive yields for pions, kaons, and protons in these data. 

MC-KLN model is a Monte Carlo version of the fac- 
torized Kharzeev-Levin-Nardi (fKLN) model [3(| ■ In the 
MC-KLN model, gluon production is obtained by nu- 
merical integration of the fc t -factorized formula [3l| at 
each transverse grid. The fluctuation of gluon distribu- 
tion due to the position of hard sources (nucleons) in the 
transverse plane is taken into account in MC-KLN. Using 



the thickness function Ta, we parametrize the saturation 
scale for a nucleus A as 



gt^;^)=2GcV 2 -^^ 
1.53 1m 



0.01 



(2) 



and similarly for a nucleus B. We choose A = 0.28 and 
a proportionality constant in the unintegrated gluon dis- 
tribution in the A: t -factorized formula to reproduce cen- 
trality dependence of px spectra for pions, kaons, and 
protons as above. 

Using the same parameter set as above, we calculate 
initial entropy distribution in U+U collisions by chang- 
ing the nuclear density from gold to uranium. To take 
account of the prolate deformation of uranium nuclei, we 
parametrize the radius parameter in the Woods-Saxon 
distribution as 

R(6, 4>) =Ro(l + hY 20 {e, <j>) + /3 4 Y 4 o(0, 0)) , (3) 

where Yi m is the spherical harmonic function, Rq = 6.86 
fm, f3 2 = 0.28 and /3 4 = 0.093 Hl|. Note that to account 
of the finite size of nucleons in the Monte Carlo approach, 
we have adjusted Rq above and the diffuseness parameter 
Sr = 0.44 to retain the nuclear density as in the original 
Woods-Saxon distribution [27[. We also take into account 
that colliding uranium nuclei are randomly oriented in 
each event. 

Figure [1] (a) shows initial eccentricity with respect 
to participant plane in Au+Au and U+U collisions at 
y/sNN = 200 GeV as a function of the number of partic- 
ipants. At each of the ten centrality bins the average ec- 
centricity and the average number of participants (AT part ) 
were calculated using both the MC-Glauber and the MC- 
KLN models. Since the eccentricity is measured in the 
participant plane, it is finite even in the very central (0- 
5%) Au+Au collisions. As previously known, the MC- 
KLN model leads to ^20-30% larger eccentricity than 
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the MC-Glauber model except in the most central events 
[30L [331 ] . In most central 5% of U+U collisions eccentric- 
ity reaches 0.146 in the MC-Glauber model and 0.148 in 
the MC-KLN model. The eccentricity is larger in U+U 
than in Au+Au collisions. Due to the deformed shape 
of uranium nucleus, this holds not only at fixed number 
of participants, but also at fixed ccntrality. The differ- 
ence, however, decreases with decreasing centrality, and 
there is almost no difference in the very peripheral events 
(70-80%). 

In Fig.[T](b), v 2 in Au+Au collisions is compared with 
the v 2 in U+U collisions. Since the rule of thumb is that 
larger eccentricity leads to larger momentum anisotropy 
and v 2 , the systematics of «2 (Apart) is similar to that 
of £ P art (Ap ar t): v 2 is larger in U+U collisions than in 
Au+Au collisions, and MC-KLN initialization leads to 
larger v 2 than MC-Glauber initialization. As well, v 2 first 
increases with decreasing A part , which reflects increasing 
initial eccentricity, but once A part falls below ~ 50, v 2 
begins to decrease. This is due to the short lifetime of 
the system which does not allow the flow to fully build 
up, and to the large fraction of the lifetime spent in the 
hadronic phase where dissipative effects are large. 

Results from the MC-Glauber initialization almost re- 
produce the PHOBOS data [34[ in Au+Au collisions. 
This indicates that there is little room for QGP viscosity 
in the model calculations. On the other hand, appar- 
ent discrepancy between the results from the MC-KLN 
initialization and the PHOBOS data means that viscous 
corrections during the plasma phase are required. 

Within the color glass condensate picture, the colli- 
sion energy dependence is taken into account through 
the saturation scale, Q s . This allows us to simulate the 
Pb+Pb collisions at y/sWJj = 2.76 TcV by using the MC- 
KLN model, adjusting the collision energy parameter and 
the nuclear density parametrization, and keeping all the 
other parameters unchanged. This is a consistent way to 
study the differences between collisions at a/s/vw = 62.4 
GeV and 2.76 TeV energies, but it may be too naive, 
since the MC-KLN model does not take into account run- 
ning coupling corrections to the evolution equation [35j . 
At RHIC energies these effects are known to have only a 
small effect, but at LHC they lead to a clearly lower mul- 
tiplicity [H, • On the other hand, these effects hardly 
affect the eccentricity [36||, which allows us to study the 
effects of the uncertainty in the final particle multiplicity 
simply by adjusting the overall factor in the unintegratcd 
gluon distribution function. Our default approach is to 
use the MC-KLN model with the same factor than in the 
RHIC calculations. In 5% most central Pb+Pb collisions 
this leads to multiplicity dN^/dr] ~ 1600 at midrapidity 
(| ?] |< 1). We also reduce the factor to obtain multiplici- 
ties dN c h/dr] ~1400 (set 1), as predicted in Ref. [Hfl, and 
-1200 (set 2)0. 



2 rj is not the shear viscous coefficient but the pseudorapidity. 



Our result for the initial state eccentricity as function 
of the number of participants in Au+Au collisions at 
^s NN —200 GeV, and in Pb+Pb collisions at y/s NN 
= 2.76 TeV is shown in Fig. [5] (a). As mentioned, the 
uncertainty in the multiplicity in collisions at ^snn = 
2.76 TcV docs not affect the eccentricity and we show 
the result obtained using our default setting. For a. fixed 
Vpart j eccentricity at LHC is apparently larger than that 
at RHIC. However, this is due solely to the larger size 
of colliding nuclei. If one compares the eccentricity at a 
fixed centrality (see each point in the figure), eccentrici- 
ties are essentially the same. 

In Fig.[2](b), v 2 in Pb+Pb collisions at ^/s^n = 2.76 
TeV is shown as a function of the number of participants 
for three different multiplicities in central collisions. The 
larger the multiplicity, the larger the v 2 , but even at the 
lowest setting of multiplicity, v 2 is clearly larger than in 
the Au+Au collisions at ^/snn = 200 GeV. 

This behavior is clearly visible in Fig. [3] where we plot 
w 2 Apart as a function of the transverse charged particle 
density (1/ S)dN c h/drj at midrapidity (| 77 |< 1) for vari- 
ous collision systems and energies. First, as expected, the 
system in U+U collisions at ^/snn =200 GeV is denser 
than in Au+Au collisions at the same energy. At ini- 
tial time to = 0.6 fm/c, the maximum temperature (en- 
ergy density) in the most central 5% of U+U collisions 
is T = 367 McV (e = 33.4 GcV/fm 3 ) and T = 361 
McV (e = 31.4 GcV/fm 3 ) in the Au+Au collisions of 
the same ccntrality. This corresponds to charged par- 
ticle transverse densities of 25.4 and 24.1, respectively, 
which means that the transverse density in U+U colli- 
sions is indeed larger, but only by — 69(0. In spite of 
the differences in the colliding systems, results for vari- 
ous centralitics in U+U collisions almost trace the ones 
in Au+Au collisions, which would suggest existence of 
scaling behavior in v 2 /e pa-r t versus (1/ S)dN c h/dri. 

However, the behavior of v 2 /e part m Pb+Pb collisions 
at ^snn = 2.76 TeV is very different. In these collisions 
the system is much denser than in the collisions at RHIC 
energies. The maximum temperature at the initial time 
T = 0.6 fm/c is T = 474, 456, 436 MeV (e = 96.2, 81.7, 
67.8 GcV/fm 3 ) for dN ch /dn ~ 1600, 1400, 1200 in 5% 
most central collisions, respectively, which corresponds 
to roughly 2-2.5 times larger transverse density than in 
most central Au+Au collisions at ^/snn = 200 GeV. 
As can be seen, v 2 /e par t no longer follows the scaling 
curve seen at the top RHIC energy, but it reaches ~ 
0.26-0.3 in central collisions. This value is — 20-35% 
larger than the value at RHIC, which is often considered 
a hydrodynamical upper limit in the literature (37j . 

The reason for this breaking of the scaling is not in 
the cascade treatment of the hadronic phase. If any- 



3 With sufficient statistics, one may make more severe centrality 
cut (e.g., 0-3%) to obtain larger transverse particle density. Mul- 
tiplicity fluctuation in the centrality cut, which we do not take 
into account, could also enhance the transverse particle density. 
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FIG. 2. (Color online) (a) e pa rt as a function of N pax t in Au+Au collisions at ^snn = 200 GeV (dashed) and in Pb+Pb 
collisions at ^/snn = 2.76 TeV (solid), (b) V2 as function of number of participants AT par t in Pb+Pb collisions at ^/snn = 2.76 
TeV for three different multiplicities in 0-5% centrality: dN^/drj ~ 1600 (solid), 1400 (dotted) and 1200 (dashed). Each point 
from right to left corresponds to 0-5, 5-10, 10-15, 15-20, 20-30, 30-40, 40-50, 50-60, 60-70, and 70-80% centrality, respectively. 
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FIG. 3. V2/e pa ,rt as a function of transverse density in 
Au+Au collisions at y / ijvjv = 62.4 (dash-dotted) and 200 GeV 
(dashed), in U+U collisions at ^/snn = 200 GeV (dotted), 
and in Pb+Pb collisions at ^snn~ 2.76 TeV (band). The 
band depicting the Pb+Pb collisions spans the results ob- 
tained using the multiplicities 1200 < dN^/drj < 1600 in 5% 
most central central collisions. 



thing, dissipation should reduce v-2, so hadron cascade 
cannot be responsible for the large value of V2 /e P art seen 
here. We have also checked that at LHC, the major part 
of V2 is generated during the hydrodynamical stage, and 
the effects of hadronic cascade are less important than 
at RHIC. The ratio of V2 generated during the hydrody- 
namical evolution to the final V2, f2,fluid/^2) in collisions 
at RHIC and LHC is shown in Fig.|4] As can be seen, the 
contribution of hadronic cascade to the total V2 with de- 
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FIG. 4. The ratio of V2 generated during the hydrodynamical 
evolution to the final V2 , f2,fluid/«2, in Au+Au and Pb+Pb 
collisions at y / sjvjv = 200 GeV and 2.76 TeV, respectively. 



fault setting at LHC is less than 5% at most centralities. 
On the other hand, the contribution reaches 10-20% of 
the total V2 in Au+Au collisions at the top RHIC energy. 

To further study the collision energy dependence of 
«2 /e P art, we do the calculation using the lower RHIC 
energy ^/snn = 62.4 GeV. As seen in Fig. |31 the ra- 
tio in central collisions (0-30%) deviates from the scaling 
curve seen at ^/sWn = 200 GeV, but the amount of the 
deviation might be too small to be experimentally ob- 
servable. The collision energy independence of i>2/e P art 
is seen in (l/S)dN c \ l /drj< 10 (fm -2 ) which corresponds 
to dNch/dr]< 50 where hadronic cascading plays a major 
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role in the whole evolution. Note that the collision energy 
dependence of W2/e pa rt is consistent with the early cal- 
culations where v<z continuously increases with total pion 
multiplicity dN v / dy at midrapidity up to 3000 at a fixed 
impact parameter (b = 7 fm) Q. See also Ref. [IH for 
a previous calculation in which the bag model equation 
of state and a higher switching temperature T sw = 169 
MeV were used. 

Recently, a prediction of elliptic flow as a function of 
transverse charged particle density up to the LHC ener- 
gies was made using viscous hydrodynamics in Ref. [39| . 
To avoid the uncertainties associated with the freeze-out 
process, vi was evaluated in that paper by calculating 
the momentum anisotropy 

= fdxdyjT^-Tvy) 
Gp J dxdy(T xx + T yy ) 1 ' 

and relying on an empirical formula «2 ~ e p /2 [40j . We 
have checked the validity of this formula in our calcula- 
tions, and found that in the collisions at the LHC energy, 
the ratio is rather vije v « 2/3, not 1/2. This discrep- 
ancy is not surprising. First, it is known that the ratio 
strongly depends on the freeze-out temperature (4lj . The 
momentum anisotropy depicts the anisotropy of the col- 
lective motion, whereas reflects the anisotropy of the 
momenta of individual particles, which includes thermal 
motion, the effects due to resonance decays [42[ and to 
the shape of the source j4^|. Second, the formula was 
found to hold in ideal fluid calculations. How dissipation 
affects it cannot be calculated, but has to be observed on 
case by case basis. 

To summarize, we predicted elliptic flow parameter Vi 



in U+U collisions at ^/sKHv = 200 GeV and in Pb+Pb 
collisions at ^/snn = 2.76 TeV using a hybrid approach 
which combines ideal hydrodynamic description of the 
QGP fluid and kinetic description of the hadronic gas. 
Due to deformation of uranium, eccentricity is larger in 
U+U collisions than in Au+Au collisions. We found the 
maximum transverse particle density is ~ 6% larger in 
0-5% central U+U collisions. «2/e P art in U+U collisions 
follows the results in Au+Au collisions, which suggests 
a scaling behavior between i^/epart and (1/ S)dN c \ l /dr]. 
However, at the LHC energy, V2/£part does not follow the 
same scaling curve and reaches the maximum value of 
~ 0.26-0.30 depending on the final particle multiplicity. 
This is clearly larger than the corresponding maximum 
value at the top RHIC energy, ^/epart ~ 0.22, and the 
so-called hydrodynamic limit for vije is not the same at 
RHIC and LHC energies. 
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